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ABSTRACT 

The application of multi- wavelength selection techniques is crucial for discov- 
ering a complete and unbiased set of Active Galactic Nuclei (AGNs). Here, we 
select a sample of 72 AGN candidates in the Extended Groth Strip (EGS) using 
deep radio and mid-infrared data from the Very Large Array (VLA) and the 
Spitzer Space Telescope, and analyze their properties across other wavelengths. 
Only 30% of these sources are detected in deep 200 ks Chandra X-ray Observa- 
tory pointings. The X-ray detected sources demonstrate moderate obscuration 
with column densities of A^h ^ 10^^ cm~^. A stacked image of sources undetected 
by Chandra shows low levels of X-ray activity, suggesting they may be faint or 
obscured AGNs. Less than 40% of our sample are selected as AGNs with optical 
broad lines, mid-infrared power laws, or X-ray detections. Thus, if our candi- 
dates are indeed AGNs, the radio/mid-infrared selection criteria we use provide 
a powerful tool for identifying sources missed by other surveys. 
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INTRODUCTION 



Active Galactic Nuclei (AGNs) are sources of tremendous energy that provide clues to 
the accretion history of the Universe. Defining a complete and reliable sample of AGNs is of 
great importance in our understanding of how black holes grow over cosmic time as well as 
the physical processes behind these powerful objects. X-ray emission is an efficient means of 
detecting energeti c activity in distant galactic cores and is often used to find AGN samples 
( iMushotzkyI |2004| ) . However, this approach is less successful at finding sources hidden by 
extreme columns of gas and dust, and models of the ob served hard X- r ay ba ckground require 
larger numbers of AGNs than are currently observed. IWorsley et al.l (120051 ) found that only 
50% of the X-ray background is resolved at > 8 keV, and the missing component of the 
background is consistent in spectral shape with a population of obscured AGNs at redshifts 
z ~ 0.5 — 1.5 with A^H ~ 10^^ — 10^^ cm~^. These synthesis models predict that at 30 keV, the 
peak of the hard X-ray spectrum, the X-ray background is dominated by heavily obscured 
sour ces and obscured AGNs outnumbe r unobscured AGNs by a ratio of three or four to 



one flUeda et all l2003l : iGiUi et all 120071 : iTreister et al.l l2005l : IComastri et al.l 120011 ). Thus 



the bulk of the accretion in the Universe is taking place in obscured environments that are 
likely to be missed by X-ray selection methods, highlighting the need for a multi-wavelength 
approach to build an unbiased census of AGNs. 

Mid-infrared and radio surveys provide a way of searching for these elusive obscured 
nuclei, as these wavelengths are much less affected by extinction and can penetrate through 
the dusty atmospheres surrounding the AGN. X-ray obscured AGNs are often bright in the 
mid-IR as the dusty environment surrounding the accretion disk reprocesses incident X- 
rays from the central engine and re-radiates these photons at infrared wavelengths. Efforts 
to select AGN candidates wit 
color criteria 



1 infrared data have ma inly concentrated on using mid-IR 



Lacv et al. 



Alonso-Herrero et al.l 



2004 



2006 



Stern et al.l 120051) or power- law selection techn iques 



Barmbv et al.l 120061 ). Recently, IPaddi et al.l (120071 ) and 



Fiore et al.l (120071 ) presented evidence of finding Compton thick AGN candidates (A^h > 10^"^) 
by searching for sources with large mid-IR excess. 

Successful attempts have a lso been made to find obscured AGNs with combined analyses 
of mid-IR and radio data sets. iDonley et al.l (|2005l ) found a minimum 4:1 ratio of obscured 
(Compton thin and Compton thick) to unobscured AGNs in the Chandra D e ep Fie ld-North 
(CDF-N) by analyzing sources with radio excess. iMartfnez-Sansigre et al.l (120051 ) found a 
2—3:1 ratio using a combination of mid-IR, near-IR, and radio fiux density criteria in the 
Spitzer First Look Survey (FLS) and the Subaru XMM- Newton Deep Field (SXDF). Both 
results are within model predictions of the obscured AGN ratio. Here, we combine and 
compare the procedures introduced by these two radio/mid-IR selection techniques to find a 
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sample of 72 AGN candidates in the Extended Groth Strip (EGS) and use multi-wavelength 
data to cross-identify their active status. In particular, understanding the X-ray properties 
of the candidate AGNs is important to control any residual contamination by non-AGNs as 
well as to study their obscuration and thus put more constraints on the X-ray background. 
Throughout this paper, we adopt A cold dark matter cosmology with Hq = 70 km s~^ Mpc~^, 
Qm = 0.3, and = 0.7. 



2. OBSERVATIONS AND DATA REDUCTION 



The Extended Groth Strip (EGS) is an approximately 2° by 15' region centered on 
a = 14^17"^ and S = -|-52°30'. The region's high ecliptic and Galactic latitudes make it 
ideal for studying extragalactic sources with little Galactic contamination, attracting deep 
coverage of the field in all wavelengths from radio to X-ray. It is a unique field in its 
combination of size and depth: it is shallower but larger than the combined GOODS fields, 
and smaller but deeper than the COSMOS field at most wavelengths. Thus, the intermediate 
position of the EGS in both depth and area makes it an id eal dataset for comb ining the AGN 
selection techniques previously tested in different regimes. iDavis et al.l (120071 ) give a detailed 
description of the region and the All- Wavelength Extended Groth Strip International Survey 
(AEGIS) data sets. 

The radio data were observed at 1.4 GHz using the Very Large Array (VLA) with ^18 h r 
of data per sky position and a resulting source sensitivity limit of 50 /iJy (jivison et al.ll2007l ). 
The data have a noise level of 10 /xJy beam~^ and a FWHM of ~ 3.8". Sources were defined 
by searching signal-to-noise ratio im ages using the SAD detection algorithm; the technique 
is outlined in lBiggs &: Ivisoru (j2006[ ). The infrared prop erties were found using instruments 



aboard the Spitzer Space Telescope (IWerner et al 



Multiband Imaging Photomet er (MIPS, 



Rieke et al 



2004), with 1.2 ks observations from the 



20041 ) and 9.1 ks exposures from the 



Infrared Array Camera (IRAC. lFazio et al.ll2004l ). The Spitzer data were processed with the 
MIPS Instrument Team Data Analysis Tool and the Spitzer Science Center IRAC pipeline, 
an d sources were identified using daophot. Further details for the MIPS data are provided 
by lEngelbracht et al.l (120071 ) and for the IRAC data by Barmby et al. (2008). The 5cr limiting 
point-source flux densities of the infrared data are 83 /jJj for the 24 //m MIPS band and 
2 /iJy for the 3.6 fim IRAC band. Their FWHMs are 5.9"and 1.8"at 24 /xm and 3.6 //m, 
respectively. 

The radio and infrared catalogs were matched to each other using a 2'.'5 search radius. 
The overlapping VLA and MIPS images result in a field of view of ~ 0.35 deg^, with roughly 
630 and 9,800 sources at 1.4 GHz and 24 ^um, respectively. In this shared area, there are 
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482 sources with both 1.4 GHz and 24 /im detections; these sources form the parent sample 
for our selection. This paper uses two radio/mi d-IR selection techniques, one of which also 
requires use of the 3.6 //m data. Due to a small offset in angle between the IRAC and MIPS 
images, the IRAC field does not extend as far north as the MIPS and radio fields. Thus, 
only 325 of the sources detected in both 1.4 GHz and 24 /zm are within the IRAC field of 
view; 321 of these have 3.6 /im counterparts. (The remaining four sources are visible on the 
IRAC image but were not extracted by daophot.) 

After selecting our candidates, we analyzed their properties using 200 ks X-ray data 
from the ACIS-I instrument on the Chandra X-ray Observatory. The X-ray image is aligned 
with the IRAC field and thus does not completely cover the VLA/MIPS shared area. Data 
reduction methods are described by Nandra et al. (2007 in prep). We also utilize R band pho- 
tometry covering the entire field to a 5a limiting AB magnitude of ~26.5 from the Suprime 
camera on the Subaru telescope (S. Miyazaki 2005, private communication). Supplemental 
R band ph otometry from th e Canada-France-Hawaii Telescope (limiting AB magnitude = 



24.2 at 8a; ICoil et al.l (120041 )) was used in a few cases where the Subaru data lacked cov- 



erage. Finally, we use redshifts obtained by the DEEP2 survey, a yailable in two-thi rds of 



the VLA/MIPS shared field from the Keck DEIMOS spectrometer flDavis et al.ll2003f ). The 
median redshift for the survey is z ~ 0.7; redshift measurements for the candidate sources 
ranged from z ~ 0.9 — 1.4. Because the optical fields are quite crowded, sources with mul- 
tiple potential counterparts within 2'.'5 were checked by eye. Two optical counterparts were 
thrown out because the match could not be determined unambiguously. 



3. AGN SELECTION TECHNIQUES 



Two different samples of radio/mi d-IR AGN candidates were selected using two separate 
sets of criteria, as illustrated in Figure 1. The first method, which we refe r to as "q-selec tion" , 
relies on the radio/infrared relation for star-forming galaxies and AGNs (ICondonlll992[). Re- 



cent f indings have suggeste d tha t the relation m a y brea k down at high redshift (IKovacs et al. 



2006; Vlahakis et al. 



20071 ) , but ISeymour et al.l (120071 ) show that the correlation holds to at 



least z ~ 3. Their study compares observed 1.4 GHz with observed 450 and 850 /im flux 
densities, which at the redshifts of interest {z > 1.5) sampl e rest wavelengt h s in t he radio 
near 5 GHz and in the far infrared near 150 and 300 /im. [Applet on et al.l (120041 ) defined 
q = log(S24 /^m/Si.4 oH z) snid showed th at q = 0.84 ± 0.28 for non-K-corrected normal star- 
forming galaxies. iDonley et al.l (120051 ) suggested that sources with strong radio emission 
relative to their 24 fim flux densities are likely to be dominated by AGNs rather than star 
formation and analyzed a sample of candidate AGNs by selecting sources with g < 0. Of the 
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482 sources detected at both 1.4 GHz and 24 fim in our field, 51 have g < 0. We discarded 
three sources that had fiux errors too large to yield convincing ratios of g < 0. A fourth 
source was discarded because the radio detection had only a marginal 24 /im counterpart 
{d = 2'.'47) and there were too many possible optical counterparts to be confident in the 
match. The q -selection only requi res the use of the VLA and MIPS parent sample. Unlike 
the sample of iDonley et al.l (120051 ) . we require all our candidates to have a detection in the 
24 i^m band in order to better compare with the second selection technique of this paper. 
Because of the offset between the MIPS and IRAC images, 35% of these sources lie outside 
the IRAC field, and 20% are outside the Chandra field. 

The g-selected samp le should include few if any starburst contaminants: adopting 
a = 0.28 at 2: = from lAppleton et al.l (120041 ). a threshold of g < selects sources 3.4 
standard deviations below the aver age g value e xpect ed for starbursts. This difference may 
even be a conservative estimate, as lBoyle et al.l (120071 ) find a higher g = 1.39 ± 0.02 for star- 



burs ts using stacked data in the Spitzer Wide Field Survey (SWIRE). iDonley et al.l (120051 ) 
and ISeymour et al.l (120071 ) show that at higher redshifts, the difference in g should remain 
constant or increase when comparing LINERs, Seyfert 2s and spiral galaxies to starburst 
systems. The g values of ULIRGs will fall closer to our selection limit with increasing red- 
shift, though their values should remain above our cut of g = out to z > 2.5. At such 
high redshifts, even our faintest candidates would have radio luminosities too high to be 
consistent with starburst systems, and we do not expect our sample to be contaminated 
with such non-AGN sources. 

The second selection method, hereafter referred to as the "flux-cut" method, employs 
the 3.6 /im data set in addition to the 24 yum and 1.4 GHz parent population to identify 
an AGN-rich sample. The selection is designed to flnd obscured AGNs near 2; ~ 2 and is 
thus a good complement to the lower-z g-selection technique. The redshift calculation comes 
from assuming the 3.6 fim e mission will be dominated by the old stellar population and 
forming a 3.6 um- z relation (jMartmez-Sansigre et al.ll2005h . analogous to the K-z relation 



for radio galaxies (IWillott et al 



I200J7 The relation assumes the galaxies have an elliptical 
host galaxy, which would not be the case for potential starburst contaminants. However, 
comparing the relation with local galaxy templates revealed that the correlation (and thus, 
the redshift assumption) remai ned valid even for reddened ULIRG templates. The flux-cut 
method is outlined in detail by iMartmez-Sansigre et al.l (120051 . 120071 ). We adopted their flux 
limits: 

(i) S24 > 300 yuJy 

(m) 100 //Jy< S1.4 GHz < 2 mJy 

{in) S3.6 /,m < 45 /uJy 
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The 24 /iin selection was chosen to find a fiux density hmited sample of active galaxies; 
a fiux density of 300 /iJy at z = 2 corresponds to quasars of luminosity L > O-'^L'^^^^^ 
at rest-frame 8 fim. At these wavelengths, dust extinction is negligible except in extreme 
cases. In AGN unification schemes, a dusty torus surrounds the optically and X-ray bright 
accretion disk. Dust particles absorb and re-emit these photons, and a to rus viewed edge- 



on wo uld t herefore be optical ly/X-ray obscured but bright in the mid-IR. IWeedman et al. 



(120061) and iBrand et al.l (l2006l ) show that the AGN fraction tends to increase as 24 /xm fiux 
density increases. Accepting only higher fiux density sources should therefore improve our 
probability of finding AGNs. 

At 2 ~ 2, the MIPS criterion alone may select star-forming systems with strong PAH 
features. The minimum radio fiux density c ut is designed to eliminate the se starburst sys- 
tems. Using the lower limit of 100 /zJy fro r nlMartmez-Sansigre et al.l (120071 ) rather than the 
350 fiJj limit of iMartmez-Sansigre et al.l (120051 ) increases our sample size while still being 
high enoug h to avoid nearly all st arburst con t amina nts. Using the equations for star-forming 
rates from lYun et al.l (120011 ) and IWu et al.l (120051 ) , if we assume that the 24 /zm emission 
of 300 /iJy arises solely from star-forming systems at z ~ 2, we would expect a radio fiux 
density that is still an order of magnitude lower than 100 /iJy. The upper limit of the radio 
criterion removes extreme radio-loud sources whose potential jets may complicate the analy- 
sis; only one potential fiux-cut candidate is omitted by this upper limit (see Figure 1). From 
the parent sample of 482 sources, 184 fit these MIPS and radio fiux criteria. 

The 3.6 /im criterion weeds out unobscured AGNs and corresponds to a rest frame 
wavelength of 1.2 fim at z = 2. Dust extinction makes the near-IR emission of obscured 
AGNs much fainter than that of unobscured sources. Requiring high 24 /im fiux densities 
and low 3.6 ^m fiux densities focuses the sample on obscured types. Out of the 184 sources 
that fit the 24 /im and 1.4 GHz criteria, 115 are within the IRAC field. Of those, the 3.6 /xm 
criterion yields a final sample size of 28 fiux-cut sources. We compared the radio/24 /xm ratios 
of these sources with the values expected from local ULIRG templates and analyzed their 
SEDs and mult i- wavelength properties. The results suggested only 2-3 potential starburst 
candidates among the sources; these few interlopers should not alter our results. 

We combined both of these techniques to arrive at our final sample. Because we are 
likely probing different redshift regimes and obscuration levels between the two techniques, 
only three sources overlap between the 47 g-selected sources and the 28 fiux-cut sources. Our 
final sample thus contains 72 radio/mid-IR AGN candidates (see Table 1). 
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4. MULTI- WAVELENGTH PROPERTIES 

4.1. X-Ray Properties 

Thirty-seven of tlie 47 g-selected sources are within the Chandra field; 15 (41%) have 
X-ray detections with Poisson false probability < 10~^ (corresponding to slightly above 3a in 
Gaussian statistics) in at least one X-ray band. Of the 28 fiux-cut sources, 26 lie within the 
X-ray image, and four have X-ray counterparts (15%). (The detection rates are only 32% and 
12% for the g-selection and fiux-cut methods in the higher significance AEGIS-X catalogs 
(Nandra et al. 2007, in prep ) , whic h require a more stringent Poisson cut of < 4 x 10~^.) 



In comparison, iDonley et al.l (120051 ) employed the g-selection method in the CDF-N, where 
VLA and MIPS observations are similar in depth to ours but Chandra X-ray exposures 
exceed 1 Ms. They found an X-ray counterpart detection rate of 40% at Poisson probability 
< 3 X 10~^ (5(7), and the rate r ises to nearly 80% for probability < 2 x 10^^ {2a) detections. 



Martinez-Sansigre et al.l (120071 ) detected > 5% of their fiux-cut-selected candidates with a 
fiux limit F2~i2 kcv = 3 x 10~^^ ergs cm~^ s~^ in the shallower SXDF with XMM exposures 
of < 100 ks. They were able to increase the detection fraction to 17% by fitting SEDs to 
filter out contaminants, a procedure we do not adopt here. 

As expected, deeper surveys are required to detect X-ray counterparts to greater com- 
pletion. Although the Chandra data in the EGS comprise the third deepest Chandra survey 
to date, 70% of our sources lack counterparts at the 200 ks depth. This may be due in part 
to the relative depth of the Spitzer data versus the Chandra field (see Figure 2). In the 
EGS, over 50% of the X-ray sources in the field have MIPS counterparts, but only 5% of 
the MIPS sources have X-ray counterparts. However, some of the undetected sources have 
upper hmits that are greater than the fiuxes of the actual detected sources, and often lie on 
relatively shallow (off-axis) areas of the X-ray map. These sources would likely have been 
detected with deeper data or on-axis pointings. 

X-ray hardness ratios are often used to study AGN type, as harder spectra generally 
indicate more obscuration. All four X-ray sources in the flux-cut method were detected only 
in the soft band, whereas 11 of the 15 sources detected with the g-selection method had both 
hard and soft band detections. Two other g-selected sources were detected only in the soft 
band, one was detected only in the hard band, and the final source was barely detected in 
the full band but had neither hard nor soft band detections at our limits. Excluding the 
source with a full band only detection, the g-selected sources have a mean hardness ratio of 
-0.19 (HR = H-S / H+S, where H = 2 - 7 keV counts and S = 0.5 - 2 keV counts). 

As measured by hardness ratios, one-thir d of the X-ray detecte d g-selected sources are 



obscured according to the division chosen by ISzokoly et al.l (120041 ) , who classify obscured 
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AGNs as those having HR > —0.2 in the Chandra Deep Field South. Our obscured frac- 
tion would rise if we include g-selected sources not detected in X-ray. Assuming all X-ray 
undetected sources are obscured would yield a maximum obscuration fraction of 73%. At 
2; = 1, we derive X-ray luminosities of lO'^^ to 10'^'' ergs/s from the observed fluxes of the 
detected sources. At this redshift, an obscuration fraction of 73% would be within the limits 



of ISzokoly et al.l (120041 ). who found 45%— 75% of their sources to be obscured in the same 



X-ray luminosity range. 

Overall, the combined AGN candidates we find tend to have softer hardness ratios than 
the general sample of X-ray sources in the AEGIS-X catalog. Including the four soft flux-cut 
sources with the g-selected sources, the mean hardness ratio falls to —0.37, while the total X- 
ray population has a mean HR = —0.17. The general population of sources with both X-ray 
and radio detections are even harder with a mean HR = —0.12. However, a lack of hard band 
detection is not necessarily an i ndication that the source has a soft spectrum, especially for 



faint sources at higher redshifts (ILaird et al.ll2006l ). Also, the g-selected AGNs generally have 



fluxes comparable to the entire X-ray sample, with an average flux F2-10 kcv = 8.3 x 10~^^ 
ergs cm~^ s~^ in the hard band and Fo.5-2 kev = 2.7 x 10~^^ ergs cm~^ s~^ in the soft band. 
On the other hand, the flux-cut sources are fainter by a factor of 10 in the soft band (and 
not detected at all in the hard band). This weaker emission, coupled with the expectation 
that these flux-cut sources should be near 2; ~ 2, suggests that the sources might be too faint 
for hard band detections, and spectral fltting is necessary to determine their true spectral 
slopes. 

Two of the X-ray counterparts with soft and hard band detections have spectroscopic 
redshifts; both of these are obscured wit h column densities > 10^^ cm~^, assuming a power- 
law models with intrinsic F = 1.9 (e.g. Nandra et al. 1997 ). Nine X-ray counterparts had 



both soft and hard band detections but no optical spectral information. We estimate column 
densities for these assuming various redshifts. At z = 0, four of the nine sources would be 
obscured with A^h > 10^^ cm~^. At 2; = 1, seven sources would be obscured, and at 2; = 2, 
eight sources would show obscuration, though none Compton-thick. The ninth source, which 
does not show obscuration, is a very bright, nearby galaxy with a hardness ratio of —0.46. 
This source is likely a normal galaxy and is further discussed in Section 4.2. 



Column densities were estimated using PIMMS (Portable, Interactive Multi-Mission Simulator; iMukai 



19931 ) v3.9b to convert between X-ray source fluxes and count rates. 
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4-. 1.1. X-Ray Stacking 
For the sources w ithout X-ray detections we use a stacking analysis to search for traces 



of X-ray emission (see iGeorgakakis et aLll2008l ). Removing sources outside the Chandra field 



and sources too close to a bright, potentially contaminating neighbor yielded 39 stackable 
sources. Stacking the counts from the 39 undetected sources yielded significant detections 
(> 3cr) in both the hard and soft bands. The mean hard band flux of the stacked galaxies 
is (1.55 ± 0.50) X 10"^^ ergs cm'^ s"^ and the mean soft band flux is (3.73 ± 0.87) x 10"^"^ 
ergs cm~^ s^^. These fluxes are near the detection limit of the 2Ms CDF-N, and deeper 
observations in the EGS would likely have detected more of these stacked sources. However, 
the average fluxes indicate that as many as half of these sources would not have been detected 
at the catalog limits even in the deepest existing X-ray survey. The work of Donley et al. 
(2005, 2007) suggests that most would be weakly detectable (> 2a) in the deepest surveys, 
with a minority well below even these limits. 

The stacked source has a hardness ratio of —0.11 ± 0.20 corresponding to an effective 
X-ray spectral slope F = 1.2, which is slightly harder than the F = 1.4 observed for the 
X-ray background. Assuming a redshift z = 1, the observed hardness ratio corresponds to a 
moderately obscured column density of ~ 4 x 10^^ cm~^. 

The distribution of detected counts in the stacked signal indicates that the emission is 
not dominated by a small sample of sources, but is representative of the overall undetected 
population as a whole. These results are for the combined sample; repeating the analysis for 
the undetected flux-cut and g-selected sources separately yields roughly similar spectra and 
fluxes, though the flux-cut sample is slightly harder with HR = —0.05 (corresponding to an 
effective F = 1.1), as opposed to HR = —0.15 (corresponding to an effective F = 1.3) for the 
g-selected sample. 

The results of the stacking analysis suggest that the undetected sources may be faint 
or moderately obscured AGNs. However, the spectral slopes are too uncertain to accurately 
determine the fraction of Compton thick sources in our sample. 



4.2. Optical Properties 

Eight of the 72 candidates have associated DEEP2 spectra, six from the g-selected 
sample and two from the flux-cut sample. Two g-selected sources (33%) and one flux-cut 
source (50%) show optical signatures of AGNs, with broad emission lines and/or strong 
absorption features indicative of outflows. All three of these optically identified AGNs also 
have X-ray counterparts. One optical source fits the spectrum of a normal massive elliptical 
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galaxy. The rest show either no optical spectral features, or features resembling star-forming 
systems. If these objects are AGNs, our selection criteria find sources that mostly have no 
optical AGN signatures and may be part of the obscured population. 

Four g-selected sources have spectroscopic redshifts ranging from 0.89 to 1.41 with a 
mean of 2; ~ 1. Two sources from the fiux-cut sample have spectroscopic redshifts, both 
near z ^ 1. Assuming the 3.6 /im emission is largely dominated by the old stellar population 
rather than AGN light, we can roughly estimate the photometric re dshift of the remaining 
sources using the K—z relation for radio galaxies (IWillott et al.ll2003l ). With our assumption, 
the 3.6 /im flux densities are upper limits to those at K, so the flux densities for the IRAC- 
detected sources imply redshifts z > 0.5. If there is an AGN contribution at 3.6 /xm (e.g., 
the power law sources discussed below), then the estimated redshifts become lower limits. If 
we assume similar host galaxies for the flux-cut sample, the se lection limit S3.6 < 45/iJy 
yields expected redshifts z > 1. iMartmez-Sansigre et al.l (120051 ) adapt the K — z relation to 
define an 5*3. e ^im — z correlation which predicts 2; > 1.4 for the fiux-cut sample. 

The hard X-ray /optical R band emission ratio (X/0; see Figure 3) is a. diagnostic widely 



used to separate AGNs from starburst galaxies (IMaccacaro et al.l Il988l : [Alexander et al. 



2OOII ) . Normal galaxies tend to have very weak hard X-ray emission and thus very small X/0 
ratios. However, AGNs are able to fuel energetic X-rays and are expected to have X-ray to 
optical ratios of ~ 0.1 — 10. An X/0 ratio of 0.01 — 0.1 is indicative of lower luminosity AGNs 
or galaxies dominated by star-formation. Of the 12 candidate sources with hard X-ray and 
R band detections in our sample, 10 have X/0 > 0.1, and 8 have X/0 > 1. Two sources in 
our sample have extreme X-ray to optical emission ratios > 10. In the general population of 
sources with X-ray and optical detections, we found that only 55% of sources had X/0 > 1. 
However, our sample sizes are too small to draw statistically significant conclusions between 
our candidates and the general population. Figure 3 shows the expected X/0 ratios from 
the stacked signals; their positions are consistent with the sources being AGNs rather than 
starburst s. 

One g-selected source, mentioned in Section 4.1, has a very small X/0 ratio and is 
extremely bright in the optical with R = 16.1. The hard X-ray fiux is F2_io kev = 1-0 x 10~^^ 



ergs cm 



This source is the only X-ray source that does not show obscuration. It is 
not plotted in Figure 3 as it is three orders of magnitude brighter in R than the brightest 
source in the figure. The X/0 ratio of this source and the soft hardness ratio suggest that 
the X-ray emission may be d ominated by X-ray binarie s and hot gas in a normal galaxy 
rather than an active nucleus (IHornschemeier et al.l 120031 ). 
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4.3. Infrared Properties 



Figure 2 plots tlie tiard X-ray flux of ttie X-ray detected sources against the 24 fim flux 
density. Half of the detected sou rces lie in the designat ed AGN region for local galaxies 
as determined observationally by iPiccinotti et al.l (119821 ) while the other half are in the 
transition region. The stacked signal is also in the transition region, tho ugh close to the 



local starbursts. Though we have not corrected our sources for redshift, [Alexander et al. 



(120011 ) show that the observed hard X-ray to mid-IR flux ratio will minimally change with 
increasing redshift for sources with low column densities and will increase more noticeably 
for obscured AGNs. 

Many AGNs can also be described by their characteristic negative power-law spectral 
energy distribu tions (SEDs) whi ch extend from mid-IR to ultraviolet wavelengths: cc u"', 
with a < —0.5 JeIvIs et al. 1994 . see also Alonso-Herrero et al. 2004 and Barmby et al. 2006). 
This emission can arise from multiple source components ne ar the dusty nucleus that broaden 



the SED rather than from a single power- law source (e.g. iRieke fc Lebofskylll98ll ). In the 
IRAC bands, stellar-dominated objects will either exhibit positive (blue) power-law behavior 
with a ~ +2 corresponding to the Rayleigh- Jeans tail of the blackbody spectrum or else not 
flt a power law at all due to their 8 /xm PAH emission or 1.6 stellar bump. A number 
of g-selected sources lie outside the IRAC observed fleld or lack detections in all four IRAC 
bands, leaving a sample of 53 out of 72 candidates for which we could study IRAC power-law 
properties. Ten of these (19%) are well flt to a power law between 3.6 pm a n d 8.0 /xm with 



a < — 0.5. E xamples of su c h flts can be found in lAlonso-Herrero et al.l (l2006l ). lBarmby et al. 
(120061 ) . and iDonley et al.l (120071 ). The detection rate is approximately the same for both 



selection methods, and is comparable t o the fraction of X-ray sources with similar power-law 
slopes in the EGS. iDonley et al.l (120071 ) show that the majority of power-law galaxies in the 
CDF-N lie at high redshifts. The low fraction of power-law detection may be due to the fact 
that most of our sample are q-selected galaxies, which are predicted to be at low-moderate 
redshifts. 

In a random sample of 53 IRAC sources, we would expect ~2 sources to exhibit negative 
IR power-law behavior based on the numbers found for the entire IRAC population. The 
Poisson probability of flnding 10 or more such power-law sources, as we have done, is only 
10~^. Of the general sample of 321 sources with radio, MIPS and IRAC detections, 23 
have IRAC power-laws. Thus, the overlap between our sample and the mid-IR power-law 
sample is substantial, and suggests that both methods are effective at flnding AGNs, but 
the populations found are distinct enough to make both methods necessary. 
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5. SUMMARY AND DISCUSSION 

We have selected 72 AGN candidates by combining and comparing two sets of radio 
and mid-IR selection criteria. The g-selection method finds 47 sources with 1.4 GHz fiux 
densities far in excess of their 24 fim fiux densities, relative to star-forming galaxies. The 
fiux-cut method selects 28 sources based on their 3.6 fim, 24 /xm, and 1.4 GHz fiux densities. 
Only three sources overlap in the two methods, since they are designed to probe different 
redshifts and obscuration levels. This selection is independent of typical selection techniques 
such as X-ray or optical identifications, and thus we can test the completeness of current 
AGN samples by finding AGNs missed at other wavelengths. 

AGNs are often defined via their mid-IR power-law behavior, optical broad lines, and 
X-ray detections. Nineteen sources (30%) in our sample have X-ray detections with Poisson 
false probabilities < 10~^ at the Chandra 200 ks depth. Column density estimates of the 
X-ray detected sources indicate that even assuming z = 0, half of our sources would be 
obscured with A^h > 10^^ cm~^. The obscuration fraction will increase with redshift. If we 
assume our sources are at z = 1, which is the mean value determined from the few sources 
with spectroscopic data, nearly all of the X-ray detected sources are obscured. Among the 
X-ray undetected sample, stacking techniques demonstrate low but significant levels of X-ray 
activity above the background emission. The stacked source has F = 1.2, which is slightly 
fiatter (harder) than the general X-ray background. The relatively high X-ray detection 
fraction and the properties of the stacked source indicate that we are finding a population of 
moderately obscured, though not Compton-thick, previously undetected AGNs. The rate of 
X-ray confirmation of AGN status is a function of the depth of the observations, and future 
work in this area will benefit from recently awarded additional Chandra pointings in the 
field. 

The radio/mi d-IR selected sample also finds more sources separately defined as AGNs 
via their optical or infrared properties than a random selection from the parent sample. 
However, our sample is also distinct from AGNs found using these other selection tech- 
nique s: less than 20% e xhibit the IRAC power-law behavior expected of luminous AGNs 
(e.g., iDonley et al.l 120071 ). and 38% show AGN signatures in their optical spectra. Table 1 
summarizes the numbers of AGNs identified by each technique separately for the fiux-cut 
and g-selected samples. Accounting for the overlap in multi-wavelength AGN identification 
for some sources, 60% of our candidates have not been identified as AGNs via optical spectra, 
IRAC power-laws, or X-ray detections. The selection criteria used here thus contribute to a 
more complete census of AGNs, and a variety of selection techniques, requiring study at all 
wavelengths, are required in order to find the least biased samples of AGNs. 



Compared to our X-ray detection rate of 30%, 



Donley et al, 



(120051) detected 80% (40%) 
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of the candidates in their (^-selected sample at 2a (Sa) in the deeper CDF-N, while iMartmez-Sansigre et al 
fl2007h detected > 5% in the shallower SXDF with a flux limit of Y2-12 kev = 3 x 10"^^ 
ergs cm~^ s~^. In general, the fiux-cut method finds a greater proportion of AGNs not 
identified via other wavelengths than the g-selected sources. This may be due to a greater 
fraction of obscured sources in the sample, as indicated by their stacked signal being slightly 
harder than the g-selected stacked signal. Also, the flux-cut sample is much fainter at X-ray 
and optical wavelengths and AGN activity may be too faint to detect in X-rays (either in- 
trinsically, or due to larger redshifts) or hidden by host galaxy light at optical wavelengths. 
However, it is also possible that the fiux-cut sample simply has a higher rate of contamination 
from starforming systems, and thus a smaller fraction of confirmed AGNs. 

Both of the selection techniques we apply find only sources that are either moder - 
ately radio-loud or occupy the bright end of the radio-quiet population (jivezic et al.ll2002l ). 
Radio-loud sources are expected to make up only ~ 10 — 20% of t he total AGN population 



flKellermann et al.lll989l : lUrrv k Padovanilll995l : llvezic et al.ll2002l ). so 80-90% of AGNs are 
undetected with these techniques. Sixty per cent of our candidates remained undetected 
using other AGN identification techniques. If we assume that under the unified model of 
AGNs, the central engine behaves the same for both radio-loud and radio-quiet sources, we 
can estimate that the same fraction of radio quiet sources are undetected by these other 
methods. This implies that half of the overall AGN population has yet to be found, and 
novel techniques must be developed in order to obtain a complete sample of AGNs. 
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Table 1. Comparison of Selection Techniques 





Number 


X-ray 


IRAC Power-law 


Optically Identified 




in Sample 


Detection 


Detection 


AGN Spectra 


Flux-Cut 


28 


4/26 


5/26 


1/2 


q- Selected 


47 


15/37 


6/28 


2/6 


Total 


72 


19/61 


10/53 


3/8 



Note. — The first column gives the number of sources selected with the radio/mid- 
IR methods outlined in this paper. The subsequent colunuis indicate the number of 
sources from our candidate list identified as AGNs via other selection techniques and 
the number of sources detected in the relevant parent sample. 
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Fig. 1. — Selection Criteria. All sources with both 24 //m and 1.4 GHz detections are plotted; 
dots represent sources with S3.6 < 45 /xJy, all other sources are depicted by crosses. AGN 
candidate sources with X-ray detections have been overplotted with squares. The flux-cut 
method selects all the dots within the dashed selection box; the g-selection method selects 
all sources above the sohd line (dots and crosses) . 
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Fig. 2. — 24 emission vs. the hard 2 — 10 keV X-ray fluxes for the sample of radio 
detected sources. Squares/large circles represent the flux-cut /g-selected sources. All other 
hard X-ray sources with 24 fim detections in the field are depicted by small dots. Filled 
symbols denote sources with X-ray detections while open symbols indicate hard X-ray upper 
limits for the non-detected sources. Filled symbols with arrows represent upper limits for 
sources detected in the X-ray, but not in the hard band. All open symbols without arrows 
are X-ray upper limits; down-facing arrows were not plotted for clarity. The two stacked 
circles/squares represent the hard band flux of the stacked signals and the mean 24 fim 
emission of the stacked sources. The horizontal line represents the flux to which 80% of the 
X-ray survey area is sensitive, and the vertical line is the 80% completeness limit for MIPS. 



Extrapolated regions for l ocal AGNs and starburst g alaxies are overplotted (jPiccinotti et al. 



19821 : iRanalh et al.ll2003l : lAlonso-Herrero et al.ll2004f ). 
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Fig. 3. — Hard X-ray (2 — 10 keV) fluxes vs. R band magnitudes. Squares/large circles rep- 
resent the flux-cut /g-selected sources. All other hard X-ray sources with R band detections 
in the field are depicted by small dots. Filled symbols denote sources with X-ray detections 
while open symbols indicate hard X-ray upper limits for the non-detected sources. Filled 
symbols with arrows represent upper limits for sources detected in the X-ray, but not in 
the hard band. Symbols with arrows pointing upper left depict candidates without firm 
detections in either hard X-rays or R band. All open symbols without arrows are X-ray 
upper hmits; left-facing arrows were not plotted for clarity. The two stacked circles/squares 
represent the hard band flux of the stacked signals and the mean R band emission of the 
stacked sources. Diag onal lines represent consta nt value X/0 ratios. AGNs are expected to 
lie above X/0 > 0.1 (IHornschemeier et al.ll2003l ) 



